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Abstract. Large-scale matter bulk flows with respect to the cosmic microwave 
background have very recently been detected on scales <100/i -1 Mpc and <300/i _1 
Mpc by using two different techniques showing an excellent agreement in the motion 
direction. However, the unexpectedly large measured amplitudes are difficult to 
understand within the context of standard ACDM cosmology. In this work we show 
that the existence of such a flow could be signaling the presence of moving dark 
energy at the time when photons decoupled from matter. We also comment on the 
relation between the direction of the CMB dipole and the preferred axis observed in 
the quadrupole in this scenario. 
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1. Introduction 

Standard Model of cosmology describes an isotropic and homogeneous universe on 
very large scales which contains (dark and baryonic) matter, radiation and dark 
energy. However, at smaller scales, the universe is no longer homogeneous since matter 
density perturbations tend to form structures that originated by primordial inflationary 
fluctuations. Thus, small volumes of matter would acquire peculiar velocities with 
respect to the CMB rest frame. These peculiar velocities are due to the gravitational 
potential of outer structures and depend crucially on the matter power spectrum. In 
standard ACDM, the large scale matter rest frame coincides with the CMB rest frame 
so that as we take larger and larger volumes of matter, peculiar velocities due to 
statistical fluctuations should become smaller and smaller. In fact, peculiar velocities 
are a very powerful tool in order to probe matter density fluctuations on very large 
scales and for that reason a great effort has been made in order to measure them by 
using different tracers: galaxies [T], clusters |2] or Type la supernovae [3]. Although 
on small scales (<60/i -1 Mpc) surveys based on different distance indicators seem to 
agree with predictions of ACDM cosmology, on larger scales they have seemed to be 
in conflict among them for many years, yielding peculiar velocities in the wide range 
— 1000 km/s. However, new data analysis performed in recent years [4J suggest that 
most of the surveys could also agree with each other on such large scales. 

A completely new approach to obtain peculiar velocities of large volumes of matter 
which makes use of the Kinematic Sunyaev-Zeldovich effect on the CMB photons by 
the hot gas in clusters of galaxies was proposed in [5] and carried out very recently in 
[6]. They find coherent bulk flows on scales of 300/i -1 Mpc towards / = 283° ± 14°, 
b = 11° ± 14° (in galactic coordinates). The reported amplitude for those peculiar 
velocities are in the range 600-1000 km/s, although the authors point out that, even 
though there is no doubt about the existence and direction of the flows, the obtained 
values for the amplitudes may have some systematic offset. In [6], they attribute these 
peculiar velocities to pre-inflationary super-Hubble perturbations (see [Tj). 

Also, in a very recent work [8], a calculation of peculiar velocities using some of the 
available measurements has been performed, but with a new method which allows to 
reduce the sensitivity to small scale power and makes possible to compare the results 
obtained from different surveys. In that work, they find a consistent flow of matter on 
scales of lOO/i" 1 Mpc towards / = 287° ± 9°, b = 8° ± 6° and with an amplitude of 
407 ± 81 km/s. The direction of the detected flow is in very good agreement with the 
results in [6], in spite of having used a very different method. The authors of this work 
claim that these peculiar velocities may be due to structures larger than the reached 
scale of 100h- 1 Mpc. 

In conclusion, these recent results suggest the existence of a coherent flow of matter 
with respect to the CMB rest frame on very large scales. In both works, such a flow is 
explained resorting to very large-scale matter perturbations, i.e., the observed peculiar 
velocities would be caused by the existence of some super-structure which must be 
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further away than the Great Attractor, located at 40 — 60/i _1 Mpc from us [9]. In 
any case, the presence of flows with such a large amplitude seems to be difficult to 
understand in the context of standard ACDM cosmology, which predicts much smaller 
velocities on the quoted scales. 

In this work, we propose an alternative explanation in which the observed bulk 
flow could be indicating that matter is globally moving with respect to radiation, i.e., 
that matter and radiation do not share a common rest frame on very large scales. 
In other words, the velocity of a matter bulk of size R may have two independent 
components: V R ,buik = V RtSta + V cosmic where V Rsta is the statistical r.m.s. velocity 
fluctuation generated by density inhomogeneities and V cosm i C is the cosmic velocity of 
matter with respect to the CMB because of having a different rest frame. Then, as we 
average over a large volume of matter, the first term becomes negligible whereas the 
second term remains constant and gives rise to a net cosmic flow. However, if that is the 
actual situation, then the question arises as how matter can be currently moving with 
respect to photons given the fact that they were strongly coupled in the early universe. 
In this work we propose that such velocity offset can be generated by the motion of dark 
energy with respect to the primordial plasma in the early universe. 



2. Moving dark energy 

In previous works [TOT [TT] , the possibility of having dark energy with a different large- 
scale rest frame from that of the matter-radiation plasma in the early universe has been 
studied. If dark energy is a perfect fluid which has always been decoupled from matter 
and radiation, there is no reason to expect it to have the same rest frame as the other 
of components of the universe. In such a case, the relative dark energy velocity vde 
should be considered as a free cosmological parameter, on equal footing with the density 
parameter DB or equation of state wde, and accordingly it could only be determined 
from observations. Indeed, the fact of having moving dark energy at the time when 
photons decoupled from matter allows to have relative velocities between dark matter, 
baryons and photons. This relative motion may be precisely the one detected in P, [8] . 

The model proposed in [TOT [TT] consists of a universe filled with four perfect fluids: 
baryons, radiation, dark matter and dark energy, whose velocities can differ from each 
other. In that case, we can define the velocity of the Cosmic Center of Mass (CCM), 
for small velocities, as: 

where a stands for the four aforementioned components of the universe. Now, if we 
place ourselves in the CCM rest frame, defined by the condition S — 0, the following 
relation must hold: 

^2(Pa +Pa)v a = 0. (2) 
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a 4 (p a +p a ) (S-v a ) =0 (4) 



In the early universe, radiation dominates over the rest of components so it will drag 
baryons and dark matter particles in such a way that they all will share a common 
rest frame, as expected for interacting species in thermal equilibrium. Concerning dark 
energy, since it does not interact with photons, it can move in a different way. However, 
as long as we use the CCM rest frame, photons and dark energy velocities are related 
by means of (j2J) as follows: 

A n early 

^arly = _^ Pr ^arly ^ 

o early . early tt \ ) 

Pde i" Pde 

where we have used that matter and radiation velocities are the same and that matter 
energy density is negligible compared to that of radiation. Notice that this relation 
provides a preferred axis in the universe given by the fluids motion direction. In other 
words, the universe will have axial symmetry around the axis defined by the motion of 
the different components. 

On the other hand, the momentum conservation equation for each fluid reads: 
d_ 
dt 

which can be immediately integrated assuming constant equation of state w a = — to 

pa 

give: 

S-v a =v°y w °- 1 (5) 

with the present value of the velocity. Thus, radiation (wr = |) moves with 
constant velocity with respect to the CCM, whereas matter (w = 0) velocity decays 
as a -1 . The evolution of dark energy velocity depends on the particular model 
under consideration. Therefore, according to (jSJ), matter and radiation have constant 
velocity before decoupling, but, after that, matter starts reducing its velocity as the 
universe expands whereas radiation keeps moving with the same constant velocity. As 
a consequence, the presence of moving dark energy at the time of recombination makes 
possible that matter and radiation could acquire a relative velocity after decoupling. 
This effect can only take place if dark energy is not a cosmological constant, otherwise 
there would be no momentum contribution from dark energy to (j2]). Notice that, since 
dark matter particles should decouple before baryons do (in order to be able to form the 
structures we observe today) and they both have the same equation of state, we expect 
them to have a relative motion with constant velocity. Thus, the complete picture of the 
motions would be as follows: radiation moving with constant velocity, baryons and dark 
matter moving in the same direction with constant relative velocity and, both, slowing 
down with respect to photons. Finally, dark energy would move along the opposite 
direction. 



3. CMB dipole from moving dark energy 



In the CCM rest frame, the metric becomes diagonal and can be written as that of an 
axisymmetric Bianchi I space-time. Moreover, as the velocities of the fluids are small, 
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this metric can be treated as a Robertson- Walker metric with a small perturbation. 
Thus, the usual expression for the dipole given by the Sachs- Wolfe effect at first order 
applies to this case, although we have to read the velocities as relative to the CCM. 
Hence, the dipole in an arbitrary frame for this model is given by [TO] : 

6 J^!±~n- (S-v)° (6) 

where n is a unitary vector along the direction of observation, v is the velocity of the 
observer today and Vd ec is the velocity of the emitter at decoupling time. Notice that 
according to ([6]) when several fluids with relative velocities are present, the CMB dipole 
is given by the velocities of emitter and observer with respect to the CCM frame. In the 
case in which all the fluids share a common rest frame, the velocites are referred to that 
frame which is nothing but the CMB frame. However, when the fluids velocities are 
different, the physically relevant frame for the dipole is the CCM which is different from 
the CMB frame. In this sense, if emission took place from a source at rest with respect 
to the CCM frame, that frame could be determined physically as the frame attached to 
an observer who measures a vanishing dipole. 

If we assume that the intrinsic dipole fluctuation at the last scattering surface is 
negligible we can take Vd ec — v^ c . On the other hand, if the observer today is at rest 
with respect to matter we have Vq ~ iP M . Thus, if we now refer the velocities to the 
CCM frame, we finally get: 



ST 



&± c n ■ (4 - i&) (7) 
J-o 

where we have used that the velocity of radiation with respect to the CCM is constant. 
Then, according to ([7]), we conclude that the contribution to the dipole temperature 
fluctuation of the CMB is due to the relative motion of matter with respect to radiation. 
This is precisely the kind of flow detected in [6l|8] so the direction reported in both papers 
gives directly the direction of motion of the fluids. 

Notice that, although calculated in the CCM frame, (J7|) is valid for any frame since 
it is expressed as the difference of two velocities evaluated at the same time. According 
to (jHJ), today the velocity of matter with respect to the CCM is expected to be much 
smaller than that of radiation and the dipole can also be written as: 

(8) 

so that the cosmological dipole could be alternatively interpreted as due to the relative 
motion of radiation with respect to the CCM. 

Using again expression (J2j), it is possible to relate the amplitude of the dipole (v^) 
to the present value of the dark energy velocity with respect to the CCM frame: 

o A 2 &DM \ , Q x 

VDE " (1 + <e)^db U R 1 + z dec + l + zj {J) 
where z* is the decoupling redshift of dark matter and radiation and w° DE is the present 
value of the dark energy equation of state. We can estimate typical values of the present 
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dark energy velocity. Thus, from the measured bulk flows we can take ~ 500 km/s 
and assuming Qde — 0.7, > 10 5 and w° DE ~ —0.97, we get v° DE ~ 1 km/s. 

4. Effects on the quadrupole 

The relative motion of the different components in the universe also affects higher 
multipoles of the CMB. In general, the contribution to the £th multipole will be ~ \v\ 
so the effect decreases rapidly as we increase I because of the smallness of the velocities. 
In particular, for quadrupolar temperature fluctuation it is possible to obtain [TT] : 

~h ( ViVj - \8ij ) nV, (10) 



T 2 V 3 

where h is a function which measures the variation of the degree of anisotropy from 
the last scattering surface until present and depends on the particular dark energy 
model under consideration and Vi is a unitary vector pointing along the direction of the 
motion of the fluids. Moreover, if such a direction is given by (6, </>), the components of 
the quadrupole are: 



a. 



20 "6^5 



fc[l+3cos20], 



'21 



[a 



A 



x )* = -J—he-^sm26, 



4 2 = {aUT = ^he-^^e. (11) 
Then, the quadrupole generated by the motion is: 

Q A , itlaiJ^k. (12) 

\ m=— 2 

On the other hand, we still have to add the isotropic quadrupole Qi coming from the 
inflationary fluctuations, whose components can be written as: 

4o = M^Qi, 



4i = -(alir=\l^e^Q 1 



42 = ^)* = ^-^^!, (13) 

where the phase factors on can be considered as stochastic variables. As the anisotropy 
is expected to be small, the total effect is given by the linear superposition of both 
contributions [12], i.e., aj m = af m + a\ m . Now, in a frame with the z-axis parallel to Vi 
we have 9 = so that the only no n- vanishing anisotropic component is a^ . That way, 
the only modified component of the total quadrupole, in that frame, is 0,20 • However, 
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in a general frame with the velocities along the direction 
becomes [TBI: 



Qt — Qa + Q 2 i ~ Z/QaQi 



where 



./ 



1 



AV5 



the total quadrupole 
(14) 

(15) 



2^6 f- 

+ 2 cos 9 cos(0 + a^)] sin 9 — (1 + 3 008(26*)) cos «i] . 

Assuming that the anisotropic quadrupole is a new contribution to add on top of the 
true model which provides the observed quadrupole, it has to satisfy [IT] : 

/iK 2 < (5T A ) 2 < 1861 /iK 2 68% C.L. 

0/iK 2 <(5T A ) 2 <5909/iK 2 95% C.L. (16) 

in order not to have a too large total quadrupole. These conditions yield the following 
constraints on the parameter h: 

0<h< 6.92 x 1(T 5 68% C.L. 



<h< 1.23 x 10" 



95% C.L. 



;i7) 



5. Model example: scaling dark energy 



Scaling models [14J are those with equation of state such that dark energy mimics the 
dominant component of the universe throughout most of the universe evolution. Thus, 
dark energy evolves as radiation before matter-radiation equality and as matter after 
that. However, in order to explain the accelerated expansion of the universe, dark energy 
has to exit from that regime and join into one with wde < — 1/3 at some point. Then, 
the evolution of the dark energy density is given by: 



Pde a T 



3w DE 



Qi e qCl 



-1 



a < a 



(:q 



Pde — \ Pde a T a 



Pde a- 3( - WDE+ ^ 



a eq < a < ax 



(18) 



a > ax 

where as commented before, a<r is the scale factor when dark energy leaves the scaling 
regime and p° DE is the present value of the dark energy density. 

In the evolution of dark energy velocity, we have to take into account the momentum 
conservation equation given by (j4j). This equation implies that the dark energy velocity 
must be discontinuous at the transition points since the equation of state jumps at those 
times whereas the quantity q^(1 ~\~w de) Pde^de is constant, being p D E continuous. With 
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this in mind, we get the following evolution for dark energy velocity in the CCM frame: 



' -/early 



VDE 



V 



DE 



4 — \->early 

gdeijO, "^DE 



a < a. 



eq 



a eq < a < ax 



(19) 



4a eo a T 3w DE ~l^farly 



3(l+w DE ) 



DE 



a > ax 
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Figure 1: Velocities evolution (CCM frame) in a scaling model with vr = 1.67 x 10 _3 c and e = 0.1. 
The continuous line (blue) is for dark energy, dashed-dotted (red) for radiation, dotted (cyan) for dark 
matter and dashed (green) for baryonic matter. Notice that in this plot dark matter is assumed to 
decouple at z* ~ 10 5 . [IT] 



In Fig. 1, we see the typical behavior of fluids velocities. In these models, 
the quadrupole can be approximated by (see [11]) Qa — 0.44 e (v^^) 2 where v^e V 
relates to radiation velocity through ([3]) as follows: vr = ev e ^ y being e the initial 
fraction of dark energy density. With this relation we can express the quadrupole as 
Q A ~ 0.44 v 2 R e~ l ~ 1.23 x 10 _6 e~ 1 or, equivalently, we have that h ~ 5.32 x KT 6 ^ 1 , 
where we have assumed that the total bulk flow is due to the motion of the fluids and 
we have taken the velocity of radiation to be vr ~ 500 km/s = 1.67 x 10~ 3 c in the 
CCM rest frame. Therefore, constraints ( 1T71) read: 

7.69 x 1(T 2 < e< 0.2 68% C.L. 



4.32 x 10 



-2 



95% C.L. 



(20) 



where the upper limit comes from primordial nucleosynthesis, which imposes the amount 
of dark energy density at that time to be less than about 20% [15] so that e< 0.2. 
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To summarize, in this kind of models it is possible to explain the presence of 
the matter bulk flow from the dark energy motion in a compatible way with the 
measurements of the CMB quadrupole. 

6. Discussion and conclusions 

Finally, we would like to comment on the fact that the total quadrupole has a preferred 
axis which happens to coincide with the direction of the velocities and, as a consequence, 
with that of the dipole. Therefore, a moving dark energy model could also shed some 
light on the so-called axis of evil problem [16J. Although this anomaly usually refers to 
the observed alignment of the I = 2 — 5 multipoles, there are also evidence that the axis 
of such alignment is correlated with the dipole direction at more than 99% CL. [17J. 
Since moving dark energy gives a common physical mechanism for both the dipole and 
quadrupole contributions, it is expected to have correlations among them. Indeed, in 
[TT] it is shown that the motion of dark energy could solve the low quadrupole anomaly 
for some models of dark energy, in particular, for scaling and null dark energy models. 
The solution of the low quadrupole problem arises because the relative motion of the 
fluids generates a certain degree of anisotropy which is seen by the photons coming 
from the last scattering surface and acquire a quadrupolar anisotropy. In a frame with 
the z-axis pointing along the fluids motions, the power of the quadrupole given by the 
generated anisotropy is zonal, i.e., it is concentrated in the m = component so that 
it gives rise to a cylindrical contribution. However, we still have to add the standard 
isotropic fluctuation generated during inflation, whose components are all comparable. 
Then, if we add linearly both contributions in such a way that the resulting quadrupole 
is lower than the inflationary one, the suppression has to take place for the m = 
component and, therefore, the total quadrupole will be non-cylindrical. 

Thus, if we want to explain the low quadrupole with this new contribution plus the 
standard inflation contribution, we need [TT] : 

54 /jK 2 < (5T A ) 2 < 3857 /iK 2 68% CL. 

0/iK 2 < (5T A ) 2 <9256/iK 2 95% CL. (21) 
which leads to the following constraints on h: 

1.18 x 1(T 5 < h < 9.96 x 1(T 5 68% CL. 

< h< 1.54 x 10~ 4 95% CL. (22) 

For scaling models, we have that the low quadrupole can be explained if the 
following constraints hold: 

5.34 x 1(T 2 < e <0.2 68% CL. 

3.45 x 1(T 2 < e <0.2 95% CL. (23) 

As the initial conditions are supposed to be random, we can compute the likelihood 
of explaining the low quadrupole with the anisotropic contribution as a function of 
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Figure 2: Probability of having a quadrupole within la interval of the measurements as a function of 
the parameter h, which contains the information about the dark energy model. 

the model-dependent parameter h. To that end, we perform simulations of the total 
quadrupole by choosing the phases and the direction (9, (ft) randomly. In Fig. 2 we 
can see that the maximum of the likelihood happens for h ~ 3.5 x 10 -5 for which the 
chance of having a quadrupole as low as the observed one is ~ 7%. 

In conclusion, we have shown that the detected large scale matter bulk flow could 
be signalling the presence of a dark energy flow which would be responsible for the 
present relative motion of matter and radiation. This kind of explanation could also 
shed some light on the problem of the low quadrupole and the alignment of dipole and 
quadrupole axes. 
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